The random growth of interfaces as a subordinated process 
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We study the random growth of surfaces from within the perspective of a single column, namely, 
the fluctuation of the column height around the mean value, y(t) = h(t) — (h(t)), which is depicted 
as being subordinated to a standard fluctuation-dissipation process with friction 7. We argue that 
the main properties of Kardar-Parisi-Zhang theory, in one dimension, are derived by identifying the 
distribution of return times to y(0) = 0, which is a truncated inverse power law, with the distribution 
of subordination times. The agreement of the theoretical prediction with the numerical treatment 
of the 1 + 1 dimensional model of ballistic deposition is remarkably good, in spite of the finite size 
effects affecting this model. 
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The random growth of surfaces is a subject of increas- 
ing interest: The number of citations of the pioneer paper 
where the Kardar-Parisi-Zhang (KPZ) equation was 
originally proposed, at the moment of writing this paper 
is of 829 in the journals of American Physical Society 
alone. The subject is discussed in excellent review pa- 
pers and books 0, Q . The interest for this field is not 
limited to the nanotechnology applications (see Refs. Q 
and for recent examples). A simple model such as 
the Ballistic Deposition (BD) [3| is an example of self- 
organization: As pointed out by Family 0, a growing 
surface spontaneously evolving into a steady state with 
universal fractal properties is similar to the mechanism of 
self-organized criticality [8( . The columns of the material 
growing due to the deposition of particles can be thought 
of as the individuals of a society. The joint action of 
the randomness driving the particle deposition and the 
interaction among columns results in the emergence of 
anomalous scaling coefficients, which can be interpreted 
as the signature of cooperation. However, only little at- 
tention has been devoted so far to studying the dynamics 
of the single individuals of this society, namely, the sin- 
gle growing columns of the sample under study. Usually 
the authors of this field of research study the correlation 
among distinct columns |IJ without paying attention to 
the dynamics of an individual. Yet, a single column is 
expected to carry information about cooperation. 

The single column perspective was recently adopted by 
Mcrikoski et al. [Tol | to study combustion fronts in paper. 
The individual property under observation is 



y(t) = h(t) - (h(t)) , 



(1) 



where h(t) denotes the height of a single column at time 
t and (h(t)) the average over the heights of the columns 
of the whole sample. The authors of Ref. 01 record the 



times at which the variable y(t) changes sign and build 
up the corresponding time series so as to create the new 
time series Tj = fcj+j — ij, namely, the set of time distances 
between two consecutive recrossings of the origin y = 0. 
The distribution density iPd{t) is shown [lfj to be an 
inverse power law with index fio, fulfilling the relation 



(3 = 2- (i D . 



(2) 



The coefficient (3 refers to the interface growth prior to 
saturation, a physical condition where the standard de- 
viation of all L columns, the interface width 



w(L,t) 
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grows as w(L,t) oc t 13 . Eq. (J2J establishes a connection 
between a single column property, /i£>, and a collective 
property, /?, thereby playing an important role for the 
perspective adopted in this paper. The theoretical foun- 
dation for this imp ortant relation is given in earlier pa- 
pers [ill IT2I ll3L IT3 and has been more recently discussed 
by Majumdar [15j . 

In this paper we prove that the KPZ condition emerges 
from the identification of iPd(t) with the distribution 
function ^(t) the essential ingredient of the subordi- 
nation theory [ilil UtI Il8| stemming from the original 
work of Montroll and Weiss I* 1 f ne subdiffusion 

case, anomalous diffusion is derived from the ordinary 
diffusion process by assuming that the time distance be- 
tween one jump and the next is determined by the inverse 
power law time distribution ipsij) with the index fi < 2. 
According to this theory 
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From the identification of [as, Eq.(4), with fir), Eq.(2), we 
obtain the anomalous scaling parameter (3 = 1/3, which 
is the KPZ prediction. 
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However, to prove that the KPZ condition is a sub- 
ordinated process, we have to show that Eq. (0 can be 
derived from the assumption that the times r, are not 
correlated, an essential property of i/>s(t). For this pur- 
pose, we assume that y is a diffusion process with scaling, 



(5) 



The number of particles located in a strip of size dy 
around y = 0, N(t), is proportional to p(0,t)dy. Thus, 
from Eq. JSJ we get 



/ s A 



(6) 



where A is a constant proportional to F(0). On the other 
hand, in the scaling regime the particles that are found 
at the origin at a given time t are only the particles that 
went back as a consequence of one or more recurrences. 
For this reason we write 



JV(t)ociJ(t) = 



(J) 



where ip n (t) is the probability that up to time t, n re- 
currences occur, the last taking place exactly at time 
t. The assumption that the times n are uncorrelated 
yields ip n (u) — (ipi(u)) n . Note that ipi(t) is the wait- 
ing time distribution function ^u(t) of which we are 
trying to assess the asymptotic properties. Note also 
that for u — > the Laplace transform of ipD(t) with 



fan - 1)T"»- 
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l /(t + ry 



is m. 



the form ipn(t) 

iJjd(u) = 1 - c u^- 1 , with c = T(2 - fi D )T' 1D - 1 . The 
parameter T < 1 is introduced to ensure the normal- 
ization condition. In the asymptotic limit of small it's, 
R(u), which is (1— cu' lD ~ 1 )/cu tJ/D ~ 1 , becomes l/cu**" -1 . 
As shown in Ref. |2l|, this is the limit for u — ► of the 
Laplace transform of a function of time, which for t — > oo 
is proportional to l/i 2- ^ . By comparison with Eq. © 
we get the important relation of Eq. (J5J) . 

The earlier remarks refer to the free diffusion part of 
the process. To derive the full picture afforded by the 
KPZ theory, free diffusion and saturation alike, we sub- 
ordinate the growth process to the ordinary Langevin 
equation 



dn 



(8) 



with the dissipation 7 fitting the condition 7 C 1 and 
T)(n) being an ordinary Gaussian white noise, defined by 



{ V (n)r)(n')) = 2Q5(\n-n'\), 



(9) 



with Q denoting the noise intensity. The dimension- 
less time I/7 3> I is essentially the saturation time of 
this ordinary fluctuation-dissipation process. Note that 
a theoretical model for the random growth of surfaces, 
called "self-similar" model by the authors of Ref. [22| . 



and proven by them to yield (3 = 1/2, is accurately de- 
scribed by the Langevin equation of Eq. © , as seen in 
Fig. ^ This means that the interaction of the column 
under study with its two nearest neighbors is properly 
taken into account by the dissipation term, i.e., the first 
term on the right hand side of Eq. JSJl . 
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FIG. 1: Theoretical standard deviation of the diffusion pro- 
cess described by Eq. JSj (solid line), compared to collapsed 
numerical results for the interface width from the "self- 
similar" model of Ref. [2^| (dots) . Q and 7 were chosen to fit 
the data points. 



At this stage we have to determine numerically the 
form of ipD{t) in the case of the BD model, which is 
known to fall in the basin of attraction of the KPZ equa- 
tion. The results are illustrated in Fig. [21 We see that 
ipo(t) is an inverse power law, with the expected index 
Hr> = 5/3, truncated at times of the order of satura- 
tion times by a faster decay. It is important to notice 
that the result of the numerical calculation illustrated in 
the insert of Fig. [3 shows that the recurrence times Tj 
are not correlated. These numerical results justify the 
identification of ifis( T ) with ipjj(r), and show that the 
subordination process must be realized with an inverse 
power law truncated at times of the order of saturation 
times by a faster decay. 

The process of subordination, yielding both anoma- 
lous free diffusion and saturation, is realized through the 
following equation 0, 0, 



dp{y, t) 



Of 



$(*-*') 7 



d , 2S f_ 

eq gy2 



dy 



y + (y 2 ) 



p(y,t')dt', 
(10) 

where $(£) is defined by means of its Laplace transform 
as 



(11) 



and (y 2 ) eq = Qh- Note that the subordination process 
has the effect of turning, in a statistical sense, time n into 
the real time t = n x l a . Thus, we have to turn the sat- 
uration time I/7 of the ordinary fluctuation-dissipation 
process into the much larger value 1/7 1 /" 3> I/7 3> 1. 
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FIG. 2: The waiting time distribution iI>(t) as a function 
of time r. The circles denote the waiting time distribution 
^d(t) of the BD model. The dashed line illustrates the theo- 
retical waiting time distribution iPs{t) used in Fig.|3]to repro- 
duce the KPZ theory, through a proper use of the BD model, 
as illustrated in the text. The insert shows the correlation 
C(t) of the waiting times Ti as a function of time r (from the 
BD model) 



At times much shorter than the saturation time 1/7 1 /", 
ips{t) is an inverse power law with index /is = /i£> = 5/3. 
This means that the Laplace transform of the memory 
kernel <I>(i) for u — > 0, for small values of u fitting though 
the condition u 3> 7 1 ^ Q , regime (i), gets the same form 
as the free memory kernel &o(t) used by the authors of 

Refs.mmm, 



3>(u) $ (u) = - it 
c 

For simplicity we set c = 

\ 1 /a 



t „.(!-<*) 



(12) 



1 , namely, T 



(l/r(l — a)) . Furthermore, we derive ips{ T ) from Eq. 
with the modified memory kernel 



$(t) = Mt)cM-i 1/a t), 



(13) 



thereby setting l>(u) of Eq. ((TT|l equal to $ (w + 7^ Q ). 
This yields V> s (w) = $ (ii + l 1/a )/{u + <J>oO + 7 1/q )), 
whose Laplace transform determines ips (t) . The compar- 
ison with the numerical iPd(t) requires some caution. In 
fact, the BD short-time behavior is model dependent and 
the subordination method is asymptotic in time. This re- 
quires that the comparison is done with a proper shift of 
tpsi T )- Nevertheless, after the shift, the decay of ips( T ) 
turns out to be slightly slower than that of ■0d( t ): an 
expected finite size effect 0- 

We now prove that the memory kernel $(i), making 
■0s(t) = ^-d( t ): compatibly with the limitation posed by 
the BD model as a substitute of the KPZ theory, satis- 
factorily reproduces the main KPZ properties. By using 
the Laplace transform approach we derive from Eq. (|10fl , 
with the form for $(i) given by Eq. I|13f) . the following 
analytical expressions for (y(t)) and (?/ 2 (i)): 



<y 2 (i)) = (y 2 (0)) *<»>(«) + {f) eq (l - K^(t)) . (15) 

The relaxation functions K^P (t) and K ^ (t) are defined 
through their Laplace transforms 



and 



k^(u) 



u + -f(u + 7 1 /a)i- 



u + 27(it + 7 1 /")i 



(16) 



(17) 



In regime (i) the functions K a (t) arc indistinguishable 
from the Mittag-Leffler functions that would result from 
setting <I>(t) = $o(i) Furthermore, for even shorter 
times correspon ding to u > 7, these Mittag-Leffler func- 
tions turn out 17] to be identical to stretched expo- 

(2) 

nentials. Consequently, K a (t) ~ exp(— 2^t a ). Mak- 
ing u larger and t shorter, the standard deviation 10 can 

be derived from Eq. (|15|) by the short-time Taylor sc- 

(2) 

ries expansion of the function Ka [t) . Thus we have 
w(t) cx t@, with (3 given by Eq. (J3J, in full accordance 
with the literature on the random growth of surfaces. 
Note that the choice of Eq. (|13J) annihilates the slow tail 
of the Mittag-Leffler function, thereby locating the onset 
of saturation in the time region of (1/7) 1 /". This can be 
easily shown by considering regime (ii), u <C 7 1 ^ Q , which 
makes Eq. (110(1 equivalent to the Fokker-Planck equation 
associated to Eq. (jHJ, with 7 replaced by 7 1 /". 

We have to point out that in real systems a regime 
of transition to the steady state exists, this being of vir- 
tually vanishing time duration in the ordinary case of 
Fig. ^ In the anomalous case this regime of transition 
becomes much more extended in time, thereby making it 
difficult to check in this case the theoretical predictions 
of Eqs. (|14fl and (|15|) . It is possible, however, to check the 
prediction of Eq. (|14f) by making an experiment at equi- 
librium, so as to avoid the out-of equilibrium induced 
aging effects. We let the system evolve till it reaches 
the steady state. Then, we label all the columns whose 
height is larger than the average height. We observe the 
time evolution of only the labelled columns and we make 
an average on many identically prepared samples. The 
result is shown in Fig. |3J and the accordance between 
theory and experiment is remarkable. 

Note that according to the current literature, the stan- 
dard deviation w(L.t), with L being the sample size, 
obeys the following form of scaling 



^2/3 



(18) 



(y(t)) = (y(0))K£\t) 



(14) 



This is valid for both the KPZ and the Edward- Wilkinson 
Q theory. f(x) is an unknown function fitting the con- 
ditions f(x) cx x@ if x <C 1 and f(x) = const if x 3> 1. 
The theory of this paper determines f{x) as follows. We 
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FIG. 3: Regression of the out-of-equilibrium fluctuation {y). 
The result of a 1 + 1 dimensional BD simulation for a system 
of size L—200 (points) is compared to the regression form 
described by Eq. 1141 (solid line), with 7 as the only fitting 
parameter, and a = 2/3. 



identify L with I/7, thereby setting x = t 2 ^^. Then, we 
note that the KPZ condition = 1/3 is obtained from the 
constraint that /j,d of Eq. J2J) is identical to fj,s of Eq. Q). 
In this condition the theory of this paper assigns to the 
unknown function f(x) the following expression 

f(x)= (l-Kf{x)) h . (19) 

To the best of our knowledge, the subordination perspec- 
tive is adopted here for the first time to account for the 
behavior of a single individual of the KPZ system. This 
is realized, on the other hand, with a proper choice of 
^s( r )> which is compatible with the transition to the 



Markov condition in the long-time limit. 

Due to the random choice of the columns the time 
distance between two consecutive arrivals of particles in 
the same column must be Poisson. How to convert this 
Poisson-like phenomenon into an anomalous subdiffusion, 
so as to involve the subordination method? The numer- 
ical results suggest that most of the arrival events are 
pseudo-events, which do not afford significant contribu- 
tions to the spreading of the diffusion distribution j^. 
Only some of these arrival events do, and the subordina- 
tion function ?/>s(t) is the distribution of time distances 
between two consecutive ones of them. 

Finally, we want to ex pres s some conjectures about the 
role of correlated noise |24j . These authors [24| proved 
that spatial and time correlation can make the scaling ex- 
ponent (3 of the EW model get values larger than (3 = 1/4. 
What about the subordination approach of this paper 
and the models with correlated noise, either EW or KPZ, 
in the d + 1 as well as 1 + 1-dimensional case? It is ex- 
pected that the time distances between two consecutive 
random events, and with it the index fis, are affected by 
both dimension and correlation. The anomalous coeffi- 
cient (3 is determined by (3 — (fis — l)/2 (EqQJ. Thus, in 
principle, it is possible to recover the value of (3 predicted 
by Ref. j2J], provided this is smaller than 1/2. Thus, the 
1 + 1 EW, with its j3 = 1 /4 can certainly be reproduced 
by means of our subordination approach. Of course, the 
condition (i$ — /J-d would be lost. This seems to be an 
interesting conclusion, since it suggests further research 
work to do to assess why fJ-s = I^d corresponds to the 
1 + 1 BD, which falls in the basin of attraction of the 
1 + 1 KPZ theory with uncorrelated noise. 
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